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We explore the use of the Inverse Amplitude Method for unitarization of scattering amplitudes 
to derive the existence and properties of possible new heavy states associated with perturbative 
extensions of the electroweak breaking sector of the Standard Model starting from the low energy 
effective theory. We use a toy effective theory generated by integrating out a heavy singlet scalar 
and compare the pole mass and width of the unitarized amplitudes with those of the original model. 

Our results show that the Inverse Amplitude Method reproduces correctly the singlet mass up to 
factors of 0(1-3), but its width is overestimated. 

PACS numbers: 


I. INTRODUCTION 

The discovery of a new particle 0,0 resembling the Standard Model (SM) Higgs boson marks the beginning of 
the direct study of the electroweak symmetry breaking sector (EWSB). The complete characterization of the EWSB 
requires the precise measurement of the Higgs couplings, as well as the search for new states. In this work we analyze 
what we can learn from the observation of departures from the SM predictions for the Higgs couplings in the case 
that no new state is observed. As is well known, anomalous Higgs couplings lead to rapid growth of the scattering 
amplitudes with energy, leading to partial-wave unitarity violation 0. Our goal is to verify how well unitarization 
procedures, more specifically the Inverse Amplitude Method (lAM) 01 predict the existence and properties of 
possible new states associated with perturbative extensions of the SM. 

Here we consider the simplest extension of the SM symmetry breaking system, i.e. we add a real singlet scalar 
field that is not charged under the SM gau ge group. Despite its simplicity, this extension of the SM can have an 
impact in the Higgs physics at the LHG [9l-ll2l|. as well as offering an interesting candidate for a portal to a hidden 
sector [l3l - [T^ . We assume that this singlet field is too heavy to be produced so we integrate it out to obtain the low 
energy effective theory. 

The lAM is based on dispersion relations to unitarize the perturbative partial wave amplitudes even in the presence 
of coupled channels, and it has been applied with success to describe low energy hadronic physics 00. This method 
has also been extensively used to study strongly interacting EWSB sectors and models exhibiting a heavy Higgs [ 13 - 
[ 21 ! • In this work we apply the lAM to the effective theory generated by integrating out a heavy singlet scalar and we 
compare its predictions to the original model parameters. In Sec.|TT]we derive the corresponding effective Lagrangian 
up to 0{p'^) and, after briefly reviewing the elements of the I AM relevant for our calculations in Sec. IHIl we present 
our results and draw our conclusions in Sec. El In particular we show that the lAM indicates correctly that only the 
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1 = 0 and J = 0 channel exhibits a resonance, reproducing the singlet mass up to factors of 0(1-3) even for relatively 
weak couplings. Its width, however, is systematically overestimated. 


II. EFFECTIVE LAGRANGIAN FOR A HEAVY SINGLET HIGGS PORTAL 

Our starting point is the SM scalar sector extended by a real singlet scalar field S 


£($,5) = + ^id^S)idf^S) -V{^,S) 


( 1 ) 


where <1? stands for the SM scalar doublet and 


V($, S) = -fiji |d>|^ + A |$|'‘ - ^ ^ ^ |$|2 S'2 ^ 


( 2 ) 


For simplicity, we have imposed a Z 2 symmetry to forbid linear and cubic terms in S. We concentrate on a scenario 
in which the S develops a vacuum expectation value (vev) vs- Presently, a new heavy scalar is allowed provided the 
ratio of the SM vev (u) to vs is small and the mixing between the mass eigenstates is small 221. 

As we will show below, in order to conveniently parametrize the low energy effective Lagrangian it is easier to write 
the SM Higgs doublet as 


$ = C/(x) ( v+H 

V2 


( 3 ) 


where 17 is a function of the goldstone bosons uJi{x) 


U {x) = exp 


;(x) 


( 4 ) 


and Ti are the Pauli matrices. Therefore, we can write Eq. o as 

C{H,S) = ^id^H){d>^H) + ^{d^S){df^S) + ^-^i±^Tr[{D^U){D^U)^ - - X^vv^H S 

Xsvs H^S + ^{2vH + H^)S^ + XvH^ + 

4 2 4 4 


( 5 ) 


with = 2Au^, M| = 2XsVg. We have traded the mass parameters and /i| for the vev’s using the the 
minimization conditions, fxjj = and /i| = Xsv^ + The covariant derivative of U takes the form: 

D^U{x) = df,U{x) + ^-gW^[x)TaU{x) - '^B^{x)U{x)t^ . (6) 

The two mass eigenstates Hi , and Si exhibit a doublet-singlet mixing due to the presence of the H S term in Eq. dSD 

Hi = cos OH -t-sin^^ and S*! = cos 0 S'— sin 0 77 (7) 

with the lighter state (Hi) identified with the recently discovered 125 GeV Higgs particle. The mixing angle 0 and 
masses are given by |23l| 

sin^ 0 =- ^ (8) 

4?/2 -I- _ 2;2 _|_ ^(1 — x2)2 + 42/2^ 

71/^2 y V 

^ (1 + T V(1-x2)2+4j/2) (9) 


with X = Mh/Ms and y = \mv/{2Xsvs)- 
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In this scenario the heavier scalar is unstable and decays via its mixing with the doublet or the singlet-doublet 
direct coupling in Eq. ©• For Mgi > 2TOtop the Si width is given by 


Ts, = r(S'i ^ W+W-) -b r(5'i ^ zz) + r(5'i -^tr) + r(5'i ^ HiHi) 


9^^Si ■ 2 r, 

7-hr- Sin d 


129,7: 


2[l-xw + -z^w ) VI - xw + xz + zx\ ) VI - xz + Sxfv/l - 


( 10 ) 


A2 


327rMsi 




XHi 


where Xi = and A/2 is the coefficient of the SiHl term obtained after we rotate Eq. ([5]) to the mass basis. Here, 


Si 


we are interested in the scenario where Si is heavy compared with Hi which allows us to approximate Eq. m by 


' 2567r 


2ll-xw + zXw VI - Xw + [l-xz + -x% VI - xz + 3a:tVl - Xt -b - xh^ 


( 11 ) 


where we have defined the parameter 


A^ 

A = -r^ so A 

As 2Ml 


~ sin^ 0. 


( 12 ) 


Si 


In the regime in which Si is very heavy we can integrate it out and generate a low energy effective Lagrangian. 
Since Hi is not a doublet field component, the corresponding effective Lagrangian cannot be expressed in terms of 
higher-dimension operators obtained in the linear representation of the electroweak symmetry breaking with a doublet 
scalar. As we will show below, it can, instead, be matched to an effective chiral Lagrangian with a light Higgs^. 

We integrate out the Si field to obtain the tree-level effective action using the approach of Ref. [1^ : the tree-level 
effective action is obtained by solving the equation of motion (EOM) and inserting the solution into the action. In 
order to do so we recast Eq. in the mass basis as 

£{Hi) + ^Si - Ml - R{x)] Si + SiB{x) + AC{Hi,Si) (13) 


with 


B{x) = ^TT[{D^U){Dfj_Ul]{Hi sin20 -b 2r!sin0) -b [Am sin20cos20 -b 2sin20(A5 sin^ 0 — Acos^ 0)] 

+ [—3sin20(Ancos0 -b Xsvs sin0) -b Am sin20(r! cos 9 + vs sin0) — Am(^^sin^ 0 -b fs cos^ 0)] 


R{x) = -^Ti[{D^U){D^Ul]sm^0 + ^H^i 


2Am (cos^ 20 - i sin^ 20 j -b 3sin2 20(A -b As) 


+Hi [Am('c cos^ 0 — vs sin^ 0) -b 3sin20(Ar!sin0 — Xsvs cos0) -b sin20Am(us cos0 — usin0)] . (14) 

A£ contains the non-quadratic terms HiSf, Sf and Sf. 

The linearized solution to the EOM for the field yields 


+ Ml + R{x) ' 

Replacing Si by Sic in Eq. (fT^ one obtains 


(15) 


CeffiHi) = C{Hi) + ]^B{x)Sic + AC{Hi,Sic) ■ 


(16) 


^ Alternatively if one integrates out the field 5, ignoring the corrections due to mixing, one can match the resulting Lagrangian to an 
effective expansion in terms of higher-dimension operators involving the remaining doublet field $ as is shown in the Appendix. 
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c 

a 

b 

d 

e 

Vc{h) 

1 

r, 


A^«^ 

Am(9Am-16AAs-10AmAs)u'‘ 

2AsM§^ 



16A^M|^ 

Vh 

0 

0 

4AsM|^ 

0 

0 

Ve 




0 

0 


8AsM|^ 

16AsM| 

Vv 

0 

0 

8AsM|^ 

0 

0 


TABLE I: Leading order in v/Ms^ coefficients defining the functions J- in Eq. (1191) . 


Now we expand the effective Lagrangian (1161) up to four derivatives and keep only terms up to dimension-six which 
allows us to match the resulting chiral effective Lagrangian to that of Refs. 

^eff{Hi) = + ccVc{Hi) + chVh[Hi) + + cjVj{Hi) - V{H,) (17) 

where 

VciHi) = ^[TTiD>^U){D^U)^TciHi) , 

Vein,) = [TTiD^^U){D^Uy]^Te{H,), 

and 

CiJ-i{Hi) = a + ^ 

We present in Table IIIl the lowest non-zero order in v/Mg^^ coefficients defining the functions T in Eq. (1191) . Within 
our approximation, the Hi potential is given by 

^ - ij) + (^ li;) ° ( 4 ) ^ 

It is interesting to notice that the operators generated at order by the integration of S*! modify the Higgs 
interactions with electroweak gauge-boson pairs {Vc) and quartic electroweak gauge-boson vertices {V&), as well as 
introducing a rescaling of all Higgs couplings to SM particles [Vh)- 


Vh{Hi) = -{d^Hi){d^Hi)TH{Hi ), 
V7 {Hi) = [Tl{D^^U){D^U)^]^,^''J^^{Hl), 


III. WEAK GAUGE BOSON SCATTERING AND ITS UNITARIZATION USING THE INVERSE 

AMPLITUDE METHOD 


The low energy effective Lagrangian in Eq. (II3 implies a modification of the gauge boson scattering with respect 
to the SM expectation, leading to unitarity violation at high energies. In this respect, two of the operators generated 
are most relevant for this discussion: Vc{h) and Veih)- Vcih) determines the Hi couplings to gauge boson pairs, in 
particular the term in ac, leads to a correction to the contribution of the virtual Hi exchange required for unitarity. 
Ve{h), in particular the term in ce, gives a contact four gauge boson coupling^. 

Eor example, the scattering amplitude at tree level for longitudinal gauge bosons is given by 


A {W+WE ^ ZlZl) = A {W+WE ^ 




Sae 


(5-2M^)(5-2M|) 


, ( 21 ) 


^ Notice despite that the UV theory is fully perturbative the effective low energy Lagrangian can be written as a theory more characteristic 
of strongly interacting/composite electroweak theories with a light scalar simply because these theories allow for enough freedom to 
account for the non-doublet nature of the light scalar. 

^ 'Pq(H) without the Higgs terms, corresponds to the L 5 operator in Refs. [2WJ^ or in Refs.[^ l2lj| . while ac and be correspond 
respectively to the coefficients 2a and h of, for example, Refs. [Uls^. 
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where y/s is the center-of-mass energy. As we can see, the term associated with ac grows as s at high energy, while 
the one containing qq exhibits growth with s^, and hence lead to violation of partial wave unitarity. 

The Inverse Amplitude Method (lAM) Q is an approach, based on dispersion relations, that allows for the full 
unitarization of the partial wave amplitudes. The lAM was originally developed for chiral perturbation theory for 
mesons and it was also applied to the unitarization of the one-loop weak gauge boson scattering amplitudes 

without a light Higgs resonance . Most recently lAM has been applied in the context of effective Lagrangians with 
a light Higgs , mostly with the aim of inferring information about the possible existence of heavier resonances 

associated with EWSB expected in composite models with a new strongly interacting sector. Let us briefly summarize 
this approach. 

The rigorous derivation of the lAM is valid only for one or several channels of particle pairs all with equal masses la¬ 
in order to apply the lAM to the longitudinal electroweak gauge boson scattering one has to work in the isospin 
symmetry approximation, i.e. setting —>■ 1 {Mz —>■ Mw = M). In this case if one defines the longitudinally 
polarized weak-gauge boson scattering amplitudes as: 

= A {Wl{py)Wl{p2) ^ fH£(P3)VLf(P4)) , ( 22 ) 

one can construct the projected amplitudes (T/) with isospin I as 

To = (0015100) = 3 A+"°°(pi,P2,P3,P4)+ A++++(pi,P2,P3,P4) , 

Ti = (10|5|10) = 2 A+"+“(pi,P2,P3,P4) - 2A+"™(pi,P2,P3,P4) - A++++(pi,P2,P3,P4) , (23) 

T2 = (20|5|20) = A++++(pi,p2,P3,P4). 

Isospin symmetry relates the amplitudes in Eq. (1^ as 

A+“ + ”(pi,p2,P3,P4) = A+“°°(P1,P2,P3,P4) + A+“™(pi, -P3, -p2,Pi) , 

A++++(P1,P2,P3,P4) = A+“°°(pi, -P3, -P2,P4) + A+"°°(pi, -P4, -P3,P2) ■ (24) 

Defining s = {pi + P 2 Y^ t = {pi — PzY = “ 4M^)(I — cos0), u = {pi — p^Y = — 5(5 — 4M^)(1 -|- cos0) with 0 

the scattering angle in the center of mass, we expand the isospin amplitudes in partial waves as: 

Ti = IOtt ^(2 J -b l)Pj(cos 0)tij (25) 

3 


where the Pj{x) are the Legendre polynomials. 

Let us assume that we know the isospin partial wave amplitudes perturbatively as 

tij = tfj -b 4 j -b ... , (26) 

where 4*/ 44 are the leading order (LO) and next to leading order (NLO) contributions. Then the lAM 

approximation (j( of the full amplitude is 


tij ~ t 


I AM 
IJ 


AO) 


(tflf 


1 _ 


(27) 


which, by construction, satisfies the unitarity constraint. 

In general, one has to deal with the possibility of coupled channels Q . For instance, in the case of chiral lagrangians 
applied to EWSB [H, [1^ the processes -A hh and hh -A hh also contribute to the partial wave I = J = 0. 

In this case one can group the perturbatively expanded amplitudes in a matrix form as 


Mqq — M, 


(0) 

00 


M, 


( 2 ) 


_ ( ‘-00 


(0) do) 


00 


A'HP ‘■HHfiOi 


^H.OO 
( 0 )' 


^.(2) ,(2) N 

'-00 '' 44,00 

h2) A2) 
Ah,0 '' 4444,002 


(28) 


where we have introduced the isospin partial wave amplitudes tnpo and Ihhpo corresponding to the processes 
W£ -A hh and hh -A hh respectively. In brief if we define 


A+ (pi,P2,P3,P4) = A{W^{pi)W^ (P2) ^ h{P3)h{Pi)) > 


(29) 
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whose projection in the 1 = 0 channel is 

Th,0 = \/3^'^“(P1,P2,P3,P4) 


the relevant partial wave is 


tH,00 


1 

647r 


1-1 


d(cos0)Tf/^o ■ 


(30) 


(31) 


And similarly tnH.oo is the corresponding partial wave amplitude for the hh —>■ hh channel. The unitarized amplitude 
matrix in this case is 


lyflAM _ ^ 

JWqq — JWqq 


(0) 


-'^^00 -'^^00 


-1 


M, 


(0) 


00 > 


so the unitarized amplitude for the WW —> WW channel is 


,IAM 

''00 




(0)',2 _ hO) 
00 > ''ff.OO 


H.oo'' 00 Vt-oo ) 


2 t 


(2) no) 


no) 


'h.oo 


( 0 ) _ ,( 2 ) _ 
''/J ''/J 


- 

Uh.oo 


nor 


j ( 2 ) -,2 

nsr 


‘‘HHfiO '‘H.OO 


(32) 


(33) 


which clearly reduces to Eq. (Ell) if the amplitude of the mixed channel {tH,oo) vanishes. 

Besides being a method for unitarization of the amplitudes, the combination of terms appearing in the denominator 
of the lAM amplitude allows for the possibility of having poles in the second Riemann sheet for some regions of the 
parameter space. When they are close enough to the physical region, those poles can be interpreted as resonances. An 
alternative approach to identify these resonances appearing in the unitarized amplitudes is to search for values of 
the center-of-mass energy (^Spoie) for which the real part of the denominator of the lAM amplitude vanishes, 

and then one identifies the mass of the resonance as = Spoie- Expanding the amplitude near the pole as 


t 


I AM 
ij 


(s) oc 


(s - M|) + 


one can also derive the value of the resonance width as oc Im[tjj“(s)]. 


(34) 


IV. RESULTS AND CONCLUSIONS 


Next we apply the lAM to unitarize the gauge boson scattering amplitudes obtained in the effective Lagrangian 
derived for the heavy singlet Higgs portal model Eq. dm. We will then search for poles in the corresponding unitarized 
amplitudes and reconstruct the properties of the inferred “resonance(s)”. In what follows we will focus on the lowest 
J partial wave amplitudes for each isospin channel, i.e. too, tn, and ^ 20 - 

Technically the mass and width of the “reconstructed” resonance are obtained by searching for poles in the denom¬ 
inator of in Eq. (l27l) i.e. by solving 


tf}{Ml)-Ret^Pj{Ml) = 0 


and 


1 Imtfj(M|) 

Mh d(^tfj{s)-Retf]{s)'^ 


ds 

Cfi 

II 


(35) 


In principle for JJ = 00 we should consider the coupled channels, which, as discussed in the previous section, are 
relevant to the WW -A WW scattering if tji.oo is not too small. For large s, tn^oo is proportional to s[bc — (ac/2)^] 
M and for the effective Lagrangian in Eq. dm this coefficient takes the value be — {acl‘^Y = which is 

Si 

assumed to be small in the effective Lagrangian expansion. So the inclusion of the coupled channels represents a small 
correction which, for simplicity, we neglect in the following and we search for the resonances in the IJ = 00 channel 
as in Eq. dm"*- The effect of the WW -A hh channel is, nevertheless, taken into account in the evaluation of ImtgQ^ 
(see Eq. (l37l) below). 


We have also verified that if we artificially set be = ciq in our calculations the reconstructed value of mass and width found in the 
7J = 00 channels are very similar to those obtained with the correct value. 
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In this work, we evaluate tree level amplitudes u sing FeynArts (3lj | with the anomalous Higgs interactions from the 
Lagrangian Eq. (IT71) introduced using FeynRules and take the exact isospin limit. Our results agree with the 
expressions in the literature in the corresponding limits. 

In order to organize the perturbative expansion of the tij we follow the counting in terms of powers of p that 
is characteristic of chiral Lagrangians [s^. In this expansion the tree-level contributions from the Higgs anomalous 
couplings, ac — 2, and be — 1, are counted as being part of tfj, i.e. 0{p'^), and therefore their corresponding loop 
contributions must be included in since they are To include the loop corrections in Ret i^T we follow the 

approach of Ref. and use the expressions obtained in Ref. [s^ using the equivalence theorem and given in 

the approximation of massless external particles. The divergent part of these loops cancel against the renormalization 
of some of the tree-level couplings of the 0{p'^) operators defined at some renormalization scale pR. This is the case 
of cg which then at a scale s becomes 


C6(s) ~ ceipji) 




{1-bc? 



(36) 


In our calculations we will take the renormalization scale as the mass of the heavy scalar pR = Thus when 

extrapolating the amplitudes to scales s ^ we can approximate C6(s) ~ CQ{Mg^) with CQ{Mg^) given in Eq. (fT51) 
and in Table I®. The remaining finite part of the loop from both the SM and the anomalous values of ac and be is 
included in as well. In order to estimate the uncertainty associated with the approximations used in the evaluation 
of this finite part we have performed our calculations both including and excluding the anomalous loop contribution. 
We will refer to these two calculations as C>(p^)-lloop and 0{p‘^)-tiee respectively. 

We first look for the presence of physical poles in the isospin amplitudes as a function of the relevant 

parameters of the effective Lagrangian: the coupling ratio A and the mass scale which determine the values 
of all relevant anomalous couplings entering the WW —>■ WW scattering, in particular ac, be and cg. One must 
notice that for the simplified potential in Eq. @ the condition that the electroweak breaking minimum is a global 
minimum sets an upper bound for A < 4A ~ 0.6; see Ref. for a recent analysis of the bounds with a more general 
potential. Nevertheless, in what follows, we will extend our study to larger values of A to illustrate the results in 
stronger coupled scenarios. 

We show in Eig.[l]contours of the real part of the denominator of the 0{p‘^)-l\oop functions i.e. Re(tj*j — 

t}J), for JJ = 00 (upper panels), 11 (central panels), and 20 (lower panels) in the s 0 Mr^ plane and for three 
characteristic values of A = 0.03, 0.3, and 3. Therefore, this figure illustrates that for no value of A do the functions 
Re(tj°^ — and Re(4o^ — t^) present a zero in the physical plane, while Re as a function of s always 

possesses a zero for any value of A and Mri. In other words, the effective theory after unitarization by the lAM 
method is compatible with the presence of one possible physical scalar resonance in the zero-isospin channel and none 
in any other spin-isospin channels, which is in agreement with the original full theory that has a scalar Si state in 
the physical spectrum and no other heavy states. 

Eor the sake of illustration, we also present in the upper panels of Fig. [1] the line corresponding to -v/s = Mr^ for 
comparison with the zero value contour which determines the position of the resonance s = Mr. As seen in this 
figure, the larger the value of A the closer the two lines, i.e. the reconstructed mass of the lAM resonance is closer 
to the real mass of the scalar of the full theory for stronger couplings. The results in the figure correspond to the 
0{p‘^)-l\oop calculation, but the same qualitative results hold for the C>(p^)-tree calculation. 

We further quantify this comparison in Fig. [5] where the upper panels depict the ratio of the reconstructed scalar 
pole mass Mr over the mass as a function of A and Mr^ for the 0(p^)-lloop calculation (left upper panel) and 
C>(p"^)-tree calculation (right upper panel). As seen in these panels, the masses agree within a factor 0(1 — 3), even 
for very small couplings independent of whether the approximate one-loop or tree amplitudes are included in the 
calculation. 

In order to verify that the scalar pole found can be interpreted as a physical state we also compute its width as in 
Eq. (1351) . For this we need to evaluate Im (s) which could be obtained from the 1-loop result in Ref. Q obtained 


® The same loops generate a coefficient for the operator [Tr(Def/)(Z)‘'f/)t][Tr(D^[/){D,/f/)t]. Such an operator is not generated by 
integrating out at the order given in Eq. lO. Thus we will take the corresponding renormalized coefficient to be zero in our 
calculations. 
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FIG. 1: Contours of the functions Re[tjj^(s) — fjj(s)] in the plane (VSjMsX for three characteristic values of the relevant 
coupling ratio A = A^/As and for the three isospin channels IJ = 00 (upper panels), IJ = 11 central panels, and IJ = 20 
(lower panels). 


using the equivalence theorem. However, since the 1-loop amplitude is only known approximately, we estimate the 
relevant imaginary part by perturbative application of the optical theorem 


Im4oHs) 






2pH 




(37) 


where p (ph) is the modulus of the three-momentum of the gauge bosons {Hi pairs) in the center of mass. 

We hnd that for all values of the model parameters A and the reconstructed width is positive, so the inter¬ 
pretation of the amplitude pole as a physical scalar state with a mass relatively close to the real scalar mass Ms^ 
holds. Notwithstanding, when compared with the perturbatively computed Si width in Eq. (mi we find that the 
reconstructed width is considerably larger as seen in the lower panels in Fig. [5J particularly for the more weakly 
interacting scenarios. This is somehow not unexpected. The lAM method was built to unitarize strong interaction 
amplitudes for which the resonance-dominance approximation holds and the amplitude near the pole of a resonance is 
fully determined by the resonance mass and width. However, for weakly interacting scenarios, such as that used here 
for illustration, the violation of unitarity is relatively mild and the full amplitude, even near the new state, contains 
a non-negligible “continuous” contribution from the SM. So the factorization used in Eqs. (|26l) and (l27l) with the full 
SM contribution included in the reconstructed amplitude as part of the resonance amplitude does not seem to be 
optimum. 

In summary, in this work we have explored the capability of the Inverse Amplitude Method for unitarization of 
scattering amplitudes to predict the properties of possible new heavy states associated with perturbative electroweak 
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0(p^)~1 loop 
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FIG. 2: Upper panels: Contours of the ratio of the mass of the resonance found in the too channel, Mr, to the mass of the 
integrated out scalar, Ms^, versus the relevant ratio of Yukawa couplings A = A^/As and Lower panels: Contours the 

ratio of the width of the resonance found in the too channel, Fh, to the width of the scalar, in the plane A (g) Msi- 


breaking extensions of the SM, using as a starting point the unitarity violating amplitudes of the low energy effective 
theory. We have used as a study case that of the singlet Higgs portal. First in Sec. [iTl we derived the effective 
Lagrangian obtained after integrating out the heavier scalar while leaving the lighter scalar, a mixture of the doublet 
and singlet states. We showed that in this case the effective Lagrangian can be matched to that of a chiral expansion 
which we write up to 0{p'^). With this effective Lagrangian in hand we obtained the relevant unitarity violating 
amplitudes. Working in the isospin approximation, we used the lAM method to reconstruct unitarized amplitudes 
and search for possible physical poles in these amplitudes. The results in Sec. IIVI show that only the unitarized spin 
scalar zero-isospin amplitude presents poles in the physical plane, in agreement with the full theory which has only 
one additional heavy scalar. We also find that the lAM reconstructs correctly the scalar singlet mass up to factors of 
0(1-3) even for relatively weak couplings. Nevertheless its width is systematically overestimated. 
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Appendix A: Effective Lagrangian after Integrating out the S Field 


Our starting point is the SM extended by the addition of a real singlet scalar field S as given in Eqs. © and ©■ 
Below the scale at which S acquires a vev Vg the scalar potential can be written as 




(Al) 


with fijj = — {XmVg)/2 and M| = 2Asu| . Now we rewrite this Lagrangian as £ = £($) J- A£($, S) with 



(A2) 
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where 


fj, = 6 Xgvs , As = 6 As , k = , A = XmVs ■ 


(A3) 


In order to apply the tree level integration procedure described in Sec. | 
EOM for the S field at lowest order leading to 


I for the singlet field S we must solve the 


where we have defined 


Sc 


1 

+ MI + U 


B 




(A4) 


(A5) 


Introducing Sc in Eq. (IA2I1 and keeping the terms up to order dimension-eight, one obtains the following anomalous 
interactions 


A£eff = 


A2 

2M| 

A2 

2M| 


i4>r + 


A2 


2M| 




- fc 1$ 


A^As 

12M| 


+ fc2- 


A fj,k 


A2 / A/r 

2M| VSMJ 


(A6) 


Ml 


2M| 


A2 

2Ml 




At this point it is interesting to apply the EOM of the doublet field to the last term of the equation above to better 
observe the emergence of an anomalous quartic coupling between the electroweak gauge bosons. This is possible since 
the invariance of the physical observables under the associated operator redefinitions is guaranteed as it has been 
proven that operators connected by the EOM lead to the same 5'-matrix elements [s^- The EOM for the doublet 
field reads 


- 2A$|$P + Fferm 


(A7) 


where Fferm is a function involving fermionic fields from the Yukawa operators. Moreover, using this EOM we find 

= 2 — 2A|$|"‘] + terms containing fermionic fields (A8) 

Therefore, altogether we find that for terms involving only scalar and/or gauge bosons 


U 


/$. 


/$,2 


f<S>,4 


/$.! 


/$.7 


A£eff = -AA|$r + :^0$.3 + ^0^,5 + ^ 0$.2 + ^0^..4 + ^0^.6 + 




Ml 






(A9) 


with AA = — 





and 


Oi.2 = 

0$,6 = 
0$,7 = 
Os,i = 


A" 


M/ 

2Xs 

3A^ / 

M 7, 

TmJ V 

3M| 


2A/^ A/i^ 

M| 

Xs Mq 

2A" / 
Mj 1 

A^Xs 1 )l2 
T2MJ ^ 

4A^ ( 

Mj V' 

-4/^, - k\ — 
2Mj 


/$.3 = 
f‘S>A = 

f-S>,6 = 


f 8A^\ A\^m 

/‘I’T = -Tv?|- = -4A— 


16A/i| 


— 12A^ 

XS Mq 


A k fi 

~1aT 


£ 2A^ 

Ml - Ml" 


^ Xs 


Xs 


+ 16A2 =.(64A2-9A^)^ 


Notice that in the last column we have introduced the relations in Eq. (IA3I1 and we have expanded to the lowest 
non-zero order in fjS^/Mg. 

The effective Lagrangian in Eq. (IA9|) also leads to violation of unitarity of electroweak boson scattering. For example 
the W~^W~ -A ZZ amplitude takes the form in Eq. m with the identification (at the lowest order in inverse powers 
of the heavy mass) oc = 2 — Cq = 0 • 


















